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Abstract 

The variational method and the Hamiltonian formalism of QCD are used to derive relativistic, 
momentum space integral equations for a quark-antiquark system with an arbitrary number of 
gluons present. As a first step, the resulting infinite chain of coupled equations is solved in the 
nonrelativistic limit by an approximate decoupling method. Comparison with experiment allows 
us to fix the quark mass and coupling constant, allowing for the calculation of the spectra of 
massive systems such as charmonium and bottomonium. Studying the results with and without 
the nonAbelian terms, we find that the presence of the nonAbelian factors yields better agreement 
with the experimental spectra. 



PACS number(s): 11.10.Ef,12.38.-t,14.40.-n 



I. INTRODUCTION 



Although the hadron spectrum is understood quite well in the context of the quark model, 
completely ab initio treatments have proved to be difficult to implement, particularly outside the 
realm of lattice gauge theory. In this paper we shall consider the description of quark-antiquark 
states using a variational approach within the canonical Hamiltonian formalism of QCD. This 
approach has been used with good effect for describing relativistic two and three body bound and 
quasi-bound states in various QFTs, including QED [1-6], the Wick-Cutkosky model [7], model 
QCD [8], and other models. An overview of the variational approach in QFT up to 1988 is given 
in the conference proceedings mentioned in [1]. A brief review of the description of few-particle 
bound and quasi-bound states in QFT by means of the variational approach in the Hamiltonian 
formalism is given in reference [9]. 

The present work is similar, in some respects, to the approach used by Zhang and Koniuk [8], 
except that they used a model Hamiltonian in which a scalar confining potential is put in by hand. 
In the present work we use the full QCD Hamiltonian in the Lorentz gauge as explained in Section 
II. The variational trial state for the quarkonium system is a Fock-space expansion that includes 
terms of the type \<j)o qq), |</>i q q g), \(f>2 q q g g), where the (j>p are variational coefficient functions 
and q, q and g represent quarks, antiquarks and gluons respectively. We use a Gupta-Bleuler type 
constraint to implement the Lorentz gauge condition. 

The variational principle then leads to an infinite chain of coupled integral equations for the 
functions <j>p, which are given in Section III. Such a system of equations is impossible to solve, so the 
rest of the paper deals with some approximate solutions. Thus, in Section IV we consider the limit 
of heavy quark masses, and approximate solutions of the nonrelativistic equations are discussed in 
Section V, as would be adequate for the heavy cc and bb systems. A comparison with some observed 
charmonium and bottomonium states is presented and discussed. Concluding remarks are given in 
Section VI. 

II. LAGRANGIAN, HAMILTONIAN AND CANONICAL QUANTIZATION 
Suppressing the flavor indices, the QCD Lagrangian density is [11] 

C Q cd = - X -F a ^ F« v + (Djij ^ - m ^ V* - (mo - m) ^ ^ (1) 

where 

F% = d»Al- d v A; + g s f abc ^ At, (2) 

and 

(DJi^Sijdn-igtfyAl. (3) 

The QCD coupling constant is denoted as g s , f a bc are the structure constants, the fermion mass 
m is an adjustable parameter not necessarily equivalent to the bare fermion mass m , and 

n = ^, (4) 

where the As are the Gell-Mann matrices. As usual, repeated indices are summed over, with the 
color indices i,j = 1, 2, 3 and a,b,c = 1, 8 for the gluon fields. 

Upon expansion of the terms above, the Lagrangian density may also be written as 

£-QCD = A/j + £>g + £>g V> + £3 g + £4 g (5) 
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where 



C$ = ft (i 7 M -m)ip % - (mo - m) ip l tp l (6) 
C g = ~d li Al(& l A av -ff'A a ' i ) (7) 

c gi , = 9s t^ft Jfft (8) 

£ 3s = -^ a "s./.ic4^ (9) 
1 
4 

Since the Lagrangian yields 



C ig = --glf abc f ade A^AlA d »A^ (10) 



n°(x) = = (11) 

we cannot quantize unless we use the Lorentz condition 

df, A a " = (12) 
and add a gauge fixing term [12] to C g such that 

C g ^~a il Al& l A av (13) 

We will use in the quantized theory, in place of (12), the weaker Gupta-Bleuler (GB) constraint 
[12] 

d^A a » + (x)\ft = 0. 
The Hamiltonian density is canonically derived from 

Hqcd = n fe ft + n a a K - c QCD (14) 

where 

n fe = OCqcd = T ii o s u aa = dC Q CD = -A aa - g s f abc A b n A ca (15) 
so that using (5) we can write 

7~Lqcd = Hip + 1~ig + H-ip g + H-ig + 7i/t g (16) 



where 



and 



= ft H l l di + m)ft + (m - m)ft ft, (17) 

H g = ~ K ^ av ~ \ di Al d l A a (18) 

H^ g = -g s t? j ft4 a ft, (19) 

n 3g = -g s f abc WA b Al +g s f abc d l A a » A\ Al, (20) 

K ig = -\g 2 s f abc f ade AlA e » [A b A d0 + A\ A d l ] . (21) 



We can expand the Dirac field (t = 0, color index j) and gluon field (t = 0, gluon index a) in 
the standard Fourier decomposition: 



■0- 



•0' 



'(x)-E/ 



d 3 p 



d 3 p 



m 



(2tt) 3 u; p 



m 

(2vr)3a; p 



[c? (p, s) u(p, s) P' x + # t (p, s ) v(p, s) e-* p ' x ] 



[d? (p, s) w(p, s) e p x + t (p, S ) n( P , a) e" 1 p ' x ] , 



a-"(x)=x; 



-7 V(2 vr) 3 2|k| 



^(k) a^(k)e lk - x + < t (k)e 



-« k-x 



and 



IP "(x) =j V ( - |k| e?(k) L?(k) e ik ' x - a" f (k) e~ 1 ^ 



(22) 
(23) 
(24) 

(25) 



where c 7 1, d J t and are the creation operators of a (free) fermion, antifermion and gluon 
respectively, and uj p = \/p 2 + m 2 . Furthermore, the anticommutation and commutation relations 
are 

{c l (p, s), c? t( p ', s ')} = R(p, s),dP Hp', s')} = <5, j 5 S „, 5 3 {p - p') (26) 
K(k),a^(k')]=CA<5 ab <5 A v <5 3 (k-k') 



(27) 



where Ci = C2 = C3 = 1> Co = — 1, and ah other commutation relations vanish. 

Since we are not interested in the vacuum energy, we will normal-order the Hamiltonian 



: Hqcd '■= j d 3 x : TCqcd ■ 



(28) 



and expand its component parts using equations (17) to (25) in (28). 

III. COUPLED EQUATIONS 

For a quark- antiquark system and an arbitrary number of gluons, we will use a Fock-space 
variational ansatz (rest frame P 'total = 0) that is a linear combination of a quark-antiquark state, a 
quark-antiquark-gluon state, a quark-antiquark-gluon-gluon state, ad infinitum. Without the terms 
that contain Fock states with at least two gluons, we would not be able to sample the nonAbelian 
terms of the Hamiltonian. Explicitly, 

\i>)= Yl J ' d 3 pd 3 nMP,n,s u s 2 ) c i \p,s 1 ) d i \n,s 2 ) <5 3 (p + n)|0) 

i,si,s 2 

+ Yl Yl J d 3 pd 3 nd 3 q 1 (f) 1 (p,n,s u S2,ci 1 ,X 1 )c i \p,s 1 )t z i )d j \n,S2)al\\ci 1 ) 
x S 3 (p + n + qi )|0) 

+ Y Y J d3p d * n d3qi d3q2 ^2(p,n,si,s 2 ,qi,q 2 , Ai, A 2 ) c l f (p,.si) 

i,j,h,si,S2 zi,Z2,Ai,A2 

x (€ h t? i3 + t? h t£ .) d> t( n , s 2 ) < t(qj) all *M S '(P + n + <H + <b)|0> + - 



3 r 

EE I d3 P d * n c l Hp,s 1 )d^(n,s 2 ) 

31, S2 i,j 

oo 3 3 8 . 

EE EE d 3 q 1 ---d 3 q /3 4> /3 (p,n,s 1 ,s 2 ,q,X) 61^61^. 

a n A \ _ 7 1 _ - „ J 



/3=0 Ai,...,A^ zi,...,Z/3' 

•it/ 



t, 



x K i a 1, • ■■%h +1 y < ' (qO a Z 1 • • • a Z fa) 5 (p + n + * + • • • + <tf)l°>» ( 29 ) 

where the A^'s are summed from to 3, i, j, k, I, m and the l^s from 1 to 3, and the z^s from 1 to 
8, and the {• • •} represent all (/?!) permutations of the Gell-Mann matrices. The total momentum 
is zero in the rest frame and the annihilation operators cP , d J and a\ have the property 

^|0> = cF|0)=<|0> = 0, 

where |0) is our trial vacuum state. In our ansatz, we have used the simplified notation [13] 

g = qi,...,q/3, A = Ai,...,A ( g 

and define qp as the set of all the g's excluding qp. We will use this same definition for the z^s and 
t Zk 's in what follows. The symmetry of <j>p{p, n, s\, s 2 , qi, . . . , c{p, Ai, . . . , A^) is apparent under the 
exchange of any (qfc,Afc) pair. 

The GB condition removes the redundant gluon degrees of freedom by requiring that our trial 
state \ip) obeys the condition 

^ = 3(q 1 )-4u(qi)^l^) = o 

for all qi, which forces the variational coefficients <pp to obey 

0i(p,n,si,s 2 ,qi,Ai = 3) = -0i(p, n, si, s 2 , qi, Ai = 0) 

2 (p,n, si,s 2 ,qi,q 2 , Ai, A 2 = 3) = -0 2 (p, n, si, s 2 , qi, q 2 , Ai, A 2 = 0) 

^ 2 (p,n, si,s 2 ,qi,q 2 , Ai = 3, A 2 ) = -<fo(p, n, si, s 2 , q 1; q 2 , Ai = 0, A 2 ) 

etc. 

The variational coefficients (f>o, <pi,...(pp in (29) are determined from the variational principle 

{5 V| : Hqcd -E:\iP)=0 (30) 

which, with our trial state and the GB condition, gives the infinite chain of coupled, multidimen- 
sional integral equations: 



uJp + Un + m (m - m) ( — + — J + |q fc | - E 
+ «/m J] Yl / - 



I 13 (f)p(p,n,s 1 ,s 2 ,q,X) 



A"=0 
2 



A"=0 s 



Qi(/3,n,s 2 ,s, A",k) 



Q 2 (p,p,s 1 ,s,X",k) 
5 



-i(p,n - k,si,s,g,k, A, A'' 



-i(p - k > n,s,s 2 ,g,k,A, A") 



/ Qzifi, k, n, S2, s) ~ 

+ ff ra LL 7r=\ </>/3-i(p,n + q fe ,si,s,g fc ,A fc ) 

s k=1 y/ml 

, sr^ ST* Q4(^, A;,p,si,s) _ ~ 

-9 m^L 7P=r (A/3-i(p + q/c,n,s,s 2 ,gfc, A fc; 

s fc=1 v l<i*l 



V ,A'''=0 k^j ' VlQfcl |qj| |k'| Iqfe + qj- -k'| 

x ^(p, n, Sl ,s 2 , q fc + q 3 - k', k', q fe q}, A'", A', X k A-,) fc, j, A'", A', k') = (31) 

with the condition 

/3 

p + n + Qfc = 0. (32) 
fc=i 

We have defined 

/ 9s 

9 - 



and the definitions for X, I 13 and the Q's are given in APPENDIX A. The factor X is composed of 
the gluon polarization vectors, and Q depends on the Dirac spinors as well as relativistic factors; 

is a color sum dependent on the Gell-Mann matrices. We should point out that the contribution 
from the "cubic" TCs g term of equation (20) vanishes because of the properties of the Gell-Mann 
matrices and our choice of variational ansatz (29) which is symmetrized in the Gell-Mann matrices. 
Consequently, the only surviving potential terms are of a color "Coulombic" type and a "quartic" 
contribution from the 7i/± g term (equation (21)); the nonAbelian effects are present only for (3 > 2. 
If, in our ansatz (29), we remove the {• • •} permutation condition on the Gell-Mann terms such 
that 

l ii 1-2 h h 1/3 Ip+i-f h h h h 1/3 l/3+i 

then we still retain the symmetry of <f>p under the exchange of any (q^, A&) pair, but now we obtain 
a set of equations that is not as tractable as the present case; that is, it's not obvious how to 
decouple them (which we need to do to make the problem tractable) . The 7i^ g contribution is now 
present and these equations are written explicitly in APPENDIX B for completeness. 

IV. HEAVY-MASS LIMIT 

The coupled equations of (31) are obviously very difficult to solve so we will begin with a 
simplification that will provide an approximate solution. One can start with the fixed or heavy 
mass limit, which will be adequate in dealing with very heavy quark systems such as J/ip (cc) or 
T (bb) that have a largely nonrelativistic behavior. 

V p 2 

In the lowest order of — , that is, letting u) v — > m-\ in the kinetic energy terms and letting 

m 2 m 

u! p — > m in the potential energy terms, (31) becomes 



p 2 n 2 P - 

2 m + + — + 2 (m - m) + V |q*| - E 
2 m 2 m ^-^ 

fc=i 



I 13 0/3 (p, n, q, A) 



2 ,. Q 



+ 9' £ Ca- / d 3 k 

A"=0 •* 



e°A»(k) 
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x [ Jf +1 <^ +1 (p, n - k, q, k, A, A") - 4 +1 0/3+1 (P - k, n, q, k, A, A") 



=1 \Mk\ 



J 3 1 ( z k) 0/3-1 (p,n + q fc ,g fc ,A fc ) - J£ ^z*) 0^_i (p + q fc , n, g fc , A fc ) 



-j- Yl ^ ^ x '" Y / 



d 3 k> 



v,v»=o - Vlqfcl M l k 'l Iqfc+q,- -k'| 

x ^(p, n, q fc + qj - k', k', q fc q~ , A'", A', A~ fc A~) X(/3, k,j, A'", A', k') = 



where 



for A" = 0, and 



e° A „(k) = e°(k)+e°(k) 



e° A „(k) = e A „(k) 



(33) 
(34) 
(35) 



for A" = 1,2. 

In the QED fixed-mass case, the kinetic energy terms are neglected and all the Fs and J's 
in (33) collapse to 1, thus allowing the infinite sequence of coupled equations to have an exact 
solution [13]: an ansatz, (ftp, may be found (along with a specific mass renormalization condition) 
that decouples the equations. However, in the QCD case the Gell-Mann matrices hamper such an 
exact solution, but we will follow the spirit of that approach in at least getting an approximate 
solution. Thus we choose an ansatz 



773 , / U / e \M^ 



x J3 0/3-1 (p,n + q0, qp, Xp) - jf 1 (zp) 0/3-i(p + q/3, n, qp, Xp) 

substitute it in (33), and impose a mass-renormalization condition (/3-dependent) 



(36) 



a 



jf +1 J 3 (zp +1 ) J% +1 Jl{zp + 



/3+1 7/3/ 



2(m -m) = - — 



+ 



d 3 k 
W 2 



7/3 7/3+1 1 7/3 7/3+1 
This uncouples the chain of equations (33) and results in a sequence of bound-state equations 

(P,n,9, A) 

. ,0 r rfih 

(p + k,n-k,g,A) 



(37) 



p 2 n 2 
2 m+ ^— + E 



2 m 2 m 



a 
4^ 



j/3+l r /3 



( z l3+l) "2 



J^ +1 j£(zp+l] 



a 



7/3 7/3+1 
2 



7/3 7/3+1 



P 



32 7T 2 I/ 3 



2 C ' x "' Y 

k+3 



A' A"'=0 



d 3 ^ 



v^fcl l^'l l k 'l hk + q.j -k' 



p, n, q fc + q, - k', k', q fc q,, A'", A', A fc Xj) X((3, k,j, A'", A', k'^ 







(38) 



with the correct rest-plus-kinetic energy for a two-particle system. It is understood that the diver- 
gent integral in (37) is controlled by a suitable regulator (cut-off) which does not appear, subse- 
quently, in (38). 



Equation (38), for any given /?, is a momentum space Schrodinger-like equation for the station- 
ary states of the heavy quark-antiquark system; there is an obvious color "Coulomb" term present 
plus a "quartic" contribution from the Tt^g term. On the surface, this "quartic" term also seems 

to share a Coulomb-type (7775-) behavior. An explicit -— 7 or ——7 term in the "quartic" part does 

I'm l*M l*M 

not appear, a term that would signal a confining-type potential. 

Note that we have chosen a specific representation [12]: 

e £(k)=< = (1,0,0,0) 
<(k) = (0,e r (k)), r = 1,2,3 

ei-e 2 = 0, e 3 (k) = — , 

k-e r (k) =0, r = 1,2 (39) 

with 

9 2 s 

4 7T 

to simplify the color "Coulomb" term. We'll work out the terms for X in the "quartic" term shortly. 

The sequence of equations (38) represents different approximations for the qq system with 
various numbers, (3, of "spectator" gluons. Note that nonAbelian effects begin only with (3 > 2. 
Explicitly, for j3 = 0, 1, 2, our mass-renormalization conditions, bound-state equations and momenta 
constraints are 

a 4 f d 3 k , A , 

nm*-m) = -—-J w (40a) 



2 m + — E 

m 



a 4 f d 3 k , . 
■°(P> n ) " ^2 3 y rj^2 <MP + k, n - k) = 



for (9 = 0, 



p + n = 



„ , . a 7 

2 (mo — m) = - — 

v ; 2vr 2 12 



d 3 k 



(406) 
(40c) 

(41a) 



p 2 n 2 

2m + f- + - £ 

2 m 2 m 



^(P,!!,^,^)-^^ 



!>i(p + k,n-k,qi,Ai) = 



(416) 



p + n + qi = 



(41c) 



for (3=1, and 



2 (mo — m) = — 



a 10 f d 3 k 



2vr 2 21 / Ikl 2 



/ 



(42a) 
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P 2 n 2 
2m+ -f- + - E 



2 m 2 m 
a 10 f d 3 k 



2(p,n,qi,q 2 , Ai, A 2 ) 

a 10 f d 6 k 
2^ 21 J |k|2 ^(P + k,n-k,q 1 ,q 2 ,A 1 ,A 2 ) 

3a A /• cPjV 

" 32^56 v ^_ n CA ' CA "' 7 Vlqillq2l|k'||qi+q 2 -k' 



A',A'"=0 

x0 2 (p,n,q 1 +q 2 -k',k',A ,/, ,A') 



P 2 n 2 

2 m + — + 



2 m 2 m 
d 3 k 



E 



X(2,l,2,A'",A',k')+A(2,2,l,A'",A',k') 
2 (p,n,qi,q 2 , Ai, A 2 ) 



tt 10 fft^ 
~ 2^2 21 J jkja 02(p + k,n-k,q 1 ,q 2 ,A 1 ,A 

~ 32 7T 2 56 ^ ^ A "' / 

\i \">— n 17 



2J 



v,v»=o - Iqal l k 'l |qi + 02 -k' 

x 2 (p, n, qi + q 2 - k', k', A'", X')X(2, 1, 2, A'", A', k') = 0, 

p + n + qi+q 2 =0 



(426) 
(42c) 



for (3 = 2. There are no gluon kinetic energy terms in (41b) and (42b) because of the ansatz (36). 

All three equations have a one-gluon exchange "Coulomb" term, modified by color factors 
stemming from the sums explicitly written in APPENDIX A; that is, we have the coefficients |, 
j2 and |y for (3 = 0, 1, 2 respectively. The (3 = equation has only the virtual one-gluon exchange 
and its familiar | factor, while the (3 = 1,2 bound state equations have spectator gluons present. 
However it is only for (3 = 2 (more generally (3 > 2) that there appears a nonAbelian contribution 
to the interquark potential energy, and so we expect that the (3 = 2 equation is the more realistic 
representation of quarkonium from among the three {(3 = 0,1,2) approximations. It is evident, 
though, that the equations become increasingly more complicated, and so more difficult to solve, 
as (3 increases. 

V. APPROXIMATE SOLUTIONS 
A. Suppression of NonAbelian terms 
If we ignore the nonAbelian terms in equation (38), we have 

d 3 k 



P 2 n 2 ^ 

2m+^- + - E 

2 m 2 m 



(pp(p,n,q,\) - 



a 
4^ 



70?) / 



with 



P + n + ^2 q fc = 0, 



2 ^(p + k,n-k,g,A)=0 (43) 



(44) 



fc=i 



and where (see APPENDIX A) 



0+1 tPi 



7(/3) 



Jf +1 Ji{ zp+1 ) , J 2 p+i JS{z p+1 ) 



(3+1 tPi 



= 2 



1(3 iP+i 

J (3+1 

03 + 1) 3 If 



+ 



1(3 1(3+1 



1(3 I (3+1 



+ 



1(3 I (3+1 



(45) 



Furthermore, unlike similar equations in QED, we note that is not a constant but /3-dependent; 
that is, it is determined by the number of "spectator" gluons. The first few terms are 



7(0) 


= 2 


4 

- = 2 (1.3333) 


7(1) 


= 2 


^ = 2 (0.5833) 


7(2) 


= 2 


^ = 2 (0.4762) 


7(3) 


= 2 


47 

- = 2 (0.4896) 


7 (4) 


= 2 


62 

— = 2 (0.4397) 


7(5) 


= 2 


= 2 (0.4357) 
3720 v ; 



(45a) 



which reveal a screening effect caused by the spectator gluons, with j((3) apparently tending to a 
limiting value. In the QED case, we would have 7(/3) = 2 for all (3. 

Equation (43) has exact solutions for each (3, as could be found easily in coordinate space. For 
example, the (3 = case of equation (43) is presented in equation (40b) and may be rewritten as 



— e 



m 



a 



4vr 2 



7(0) 



d 3 k 



,(p + k) = 



(46) 



where 4>o(p) represents 4>o(p, n) (recall that p + n = for (3 = 0). In coordinate space, equation 
(46) transforms to the familiar 



X/l _ 07(0) 
2 fi 2 |r| 



0o (r) = e 0o (r) 



(47) 



where 



and 



A 4 



e = £ - 2 m. 

Equation (47) has the standard (modified) solutions 

M« 2 7 2 (0) 



E — 2m 



8n 2 



(48) 



The 0o (r) of (47) correspond to the momentum space wave functions 



00 (p) OC 



p 



(p 2 + a 2 ) 



2\2+l ^n-l 



C 



y -a 2 

~ l[ p 2 + a 2) 



with 



fi a 7(0) 



2 n 



(49) 
(50) 



and the C^t^'s are the Gegenbauer functions. 
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Now, for (3 = 1 (which corresponds to a term with the q q pair plus an explicit gluon; i.e., an 



H2 -like configuration), equation (43) yields equation (41b), which we rewrite as 



p 2 n 2 

+ - e 



o f d k 

ii(p, n, Ax) - — 7 (1) J — MP + k, n - k, Ai) = 



2 m 2 m 

with 0i (p, n, Ai) = 0i (p, n, qi, Ai) (recall that p + n + qi = for this case). If we assume that 

0i(p,n,Ai) = 0i(p,n) Z(Ai) 
then the Z factors out and in coordinate space (51) is 



(51) 



2 m 



Q 7(1) 
2 |ri - r 2 | 



H(ri,r 2 ) = e 0i(ri,r 2 ) 



If we transform to new coordinates 



r = ri - r 2 , 
R = ^±^ 



and let 



0i(ri,r 2 ) = 0i(r) 0i(R) 
equation (52) reduces to the form (47), with a solution 

^« 2 7 2 (1) 



e = E-2m 



8n 2 



(52) 

(53) 
(54) 

(55) 
(56) 



Here, 0i (r) is the typical hydrogenic wavefunction, and <pi (R) is a plane wave solution. In momen- 
tum space, 

<Mr) <MR) ^ <M Pr ) 5(P R ) (57) 

where 0i(p r ) is the typical momentum space hydrogenic wavefunction (see (49)). Noting the 
transformations (53) and (54), the analogous procedure in momentum space is to let 



and 



p - n 



P R = p + n. 



(58) 



(59) 



This will be of use in the next section. 

Lastly, for j3 = 2 (a system that is composed of a q q pair plus two explicit gluons; i.e., 
an i7 2 -moleculedike configuration), equation (43) yields (equation (42b) without the nonAbelian 
terms) 

a f d 3 k 

6 2 (p,n, qi )- — 7 (2) J — 2 (p + k,n-k, qi ) =0 (60) 

with 02 (Pi n ) qi) = 02 (P) n ) qi) q2) (recall p + n + q x + q 2 = 0) and we have factored out the A 
terms as in the last case. In coordinate space, equation (60) is 



p 2 n 2 
-5— H 

2 m 2 m 



1 / 2 2 n a 7( 2 ) 



2 m 



2 ri - r 2 



6 2 (ri,r 2 ,r 3 ) = e 2 (ri, r 2 , r 3 ). 



(61) 
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Note that r3 appears as a parameter in equation (61); it has no effect on e. As for the (3 = 1 case, 
we can use the transformations (53) and (54) to get 



t = E-2 m =-H^, (62) 

where the wave functions have the same structure as (3 = 1. We expect the same result to hold for 
(3 > 2, with arbitrary functions present. That is, in general, we have 

£ = £ _ 2 ™=-^^a (63, 

(with an 'inverse Bohr radius' a = — — li^l use d in the respective momentum space wave func- 

2 n 

tions). 

Let us apply this modified Balmer formula to the heavy quark mesons such as charmonium and 
bottomonium to obtain a prediction for the low-lying bound states. Using the experimental values 
for the lowest-lying IS and 25 states of cc (i.e., J/ip(lS), J/i/j(2S); [11]), we can fix two of our 
parameters; that is, we can set our mass for the charm quark (m c ) as well as the coupling constant 
a particular to this system. Once this is done, we use our expressions for E (equations (48), (56), 
and (62)) to get the predicted values for the 35, 45, 55 and 65 states. We repeat the process for 
the bb mesons, using T(15), Y(25), and present the results in Tables I and II for charmonium and 
bottomonium respectively. 

In this approximation suppressing the nonAbelian terms, we note that the predicted values for 
charmonium and bottomonium are independent of (3, the number of "spectator" gluons present. 
The mass values are constant, as can be easily seen from equation (61) for n = 1,2, as is the 
value a 2 7 2 (/3); only a changes with (3. As expected, the predicted masses are not particularly 
close to the observed ones, though the values for bb are better than for charmonium, since this 
is a heavier, more nonrelativistic system for which the mass spectrum is influenced more by the 
Coulomb potential. One would expect better results still for ti, were the data available. 

B. NonAbelian terms present: (3 = 2 case 

We now turn to the more realistic approximation in which we have an explicit manifestation of 
the nonAbelian nature of the interaction. Recalling our transformations (58) and (59), our equation 
(42b), for (3 = 2, may be rewritten in the new coordinates as 



2 i> - 

2(p r ,PR,qi,q2,Ai, A 2 ) 



„ 2 I) 

Pr . r R 

e 

2(i 8fi 



a f d 3 k 

~4^ 7 ^ J jkf ^ 2< - Pr ~ k ' P R' c H' c l2, Ai,A 2 ) 

6 a y—^ f d 3 k' 

~ 32 vr2 56 y ^ =o CA ' CA "' J Vkil lq 2 | |k'| |qi+q 2 -k'| 

x MPr, Pr, qi + q 2 - k', k', A'", \')X(2, 1, 2, A'", A', k') = (64) 
If we make the transformation 

vliil v / |q2l02(Pr,PR,qi,q2,Ai,A 2 ) = $(p r ,PR,qi,q2) z(\ 1 ,x 2 ) (65) 
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then equation (64) becomes 



Pr , "R 



^(Pr,PR,qi,q2) Z(\ 1 ,\ 2 ) 



2 ii 8fi 

a I' d 3 k 

- — 7 (2) J — $(p r -k,P R , qi ,q 2 ) Z(X 1 ,X 2 

6 a A f _ d 3 k' 

~ 32 vr 2 56 ^ ^ A ' ^ A "' J Ik 7 ! 

A',A"'=0 J 1 1 



|qi + q2 -k'| 



x $( Pr ,P R , qi +q 2 -k',k') Z(X"',X')X(2,l,2,X"',X',k') = 



(66) 



For an approximate solution, let's multiply (66) by <E>(p r , Pr, qi, q2) and integrate over p r , Pr, qi 
such that 



J d 3 Pr d 3 P R d 3 qi 
,3 



™ 2 p> 2 
Pr r R 



— € 



$ 2 (p r ,P R , qi ,q 2 ) Z(Ai,A 2 ) 



2 n 8 n 

^ d 3 Pr d 3 P R d 3 qi $(p r ,P R , qi ,q 2 ) $(p r -k,P R , qi ,q 2 ) Z(Ai,A 2 ) 
d 3 fc' d 3 p r d 3 P R d 3 qi 



6 a 
32 vr 2 56 



^ Ca' Ca 



|k'| |qi +q 2 -k'| 



A',A"'=0 

X $(p r , P R , qi , q 2 ) <D(p r , P R , qi + q 2 - k', k') Z(X"', X')X(2, 1, 2, A'", A', k') = (67) 



As a trial function, we will use the nonAbelian solutions of the last section 

$(PrPR,qi,q2) = MPr) «5(Pr) ^(qi) 



(68) 



where c (p r ) is the classic hydrogenic wave function and ip{q.i) is an arbitrary function which we 
choose to regulate the integral in momentum space. We can treat a in the wave function as a 
variational parameter to optimize our results (see (49), (50)). Integrating over P R leads us to 



/ 



d 3 p r 



2 n 



— e 



( Pr ) Z(X 1 ,X 2 ) j d 3 Ql ^ 2 ( qi ) 



a f d 3 k 

4^ 7(2) / W ^ c(Pr) ^ c(Pr " k) Z(Al ' * 2 ' 



32 vr 2 56 2 Ca ' Ca "' / 
x j d 3 qi V-(qi) 1(2, 1,2, A'", A', k 



d 3 p r ^( Pr ) Z(A"',A') 




Ai V 2 (qi) 
d 3 fc' V(-k') 



rl\2 



(69) 



with 



i.e., 



qi + q 2 + P R = qi + q 2 = 0; 



q2 = -qi 
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which simplifies the polarization vector terms in X considerably (see APPENDIX A, equation 
(A8)). Summing over Ai, A 2 in (69), we obtain 



d 3 p r 

a 
4^2 



2 fj, 



32 vr 2 56 



(Pr) J Ai ^ 2 (qi) 

/d 6 k f 
d 3 Pr MPr) MPr ~ k) J d 3 Ql ^ ( Ql ) 

E C*i CA 2 / ^e(Pr) / " 

' v— n ^ ^ 



A',A"'=0 



J d 3 Ql ^(q 1 )X(2,l,2,A ,,, ,A , ,k') 



|k'| 2 

A"'<-+Ai, A'^A 2 = 



(70) 



Summing again to remove the Ai,A2 terms still present in the "quartic" term of (70), that is, 
working out 



E E CA 1 CA 2 X(2,l,2,A'",A',k') 

A',A"'=0 Ai,A 2 =0 
using the convention of [14] 

ei(-k) = -€i(k), e 2 (-k) = c 2 (k) 
and recalling (39), we can simplify (70) to 



(71) 



(72) 



/ 



d 3 p r 



2 n 



( Pr ) J d 3 Ql ^ 2 (qi) 



4^ 7(2) / W ^ MPr) " k) / d3<11 ^ 2(qi) 



12 • 36 a 
32 vr 2 56 



d 3 Pr <Pl(Pr) 



d 3 k> V(-k') 



x { ei(k') • ei(qi) + e 2 (k') • e 2 (qi) - Ci(qi) ■ -j-r 



ei(k') ■e 2 (q 1 ) + e 2 (k') -Ci(qi) -e 2 (qi) ■ j^Tj 



d 3 qi ^(qi) 

\ 2 

k' 



(73) 



or, 



e = 



dVf|^(Pr) -^7(2) / 0^Pr^ 



(Pr ) <£c (Pr - k) | y/ 7 d V </» 2 (p r ) 



^ 7 dfcX-kO 7 dgx g 2 V(qi)/ 7 ^^(qi) 



(74) 



We will use trial functions of the scaled hydrogenic-type (equation (49)); for example, for the 
ground state we have 

MPr) " (^W' 
14 



where a is an arbitrary scale parameter. For the trial function ^>(qi), we choose a form 

^ (a ' qi)= (^W' e_ ^' 

using the same scale a as for </> c (p r ). In the second part of (74), changing the variable by letting 

it 

a 2 a a 7(2) 9 a a 



q[ = — allows us to pull out the a. That is, 
a 



m 2 n 14 n 

where we define 

,2 



Y (75) 



y = J "7 .r, J 2 .7/r! yvM1 • (^) 



/ d?i ?i 2 V' 2 (qi) 

If we optimize (75) with respect to a we obtain the optimum value 



and so the corresponding 



ma / 10 9 Y \ 
a " = — 2T^ + 14T (77) 



m« 2 / 10 9Y\ 
£ " = -^ 2T^ + 14T • (78) 



Note that the binding energy is non-zero as n — > 00, quite unlike the Abelian (QED-like) case. 
We try a function 

* M) = WTW m 

in (78) and use the 15 and 2S levels to fix a and m. Thus we get the predicted values for the 
35, 45, 55 and 65 states for cc which are listed in Table III and for bb in Table IV. Alternatively, 
trying a function of type 

V(q'i) = e"^ (80) 

produces the results in Table V and Table VI, for cc and bb respectively. These results are for the 
[3 = 2 case only and we find that the inclusion of the nonAbelian terms improves the agreement 
with experiment substantially over the results of Tables I and II which were derived using only 
the modified Balmer formula (nonAbelian terms suppressed). Once again, the bottomonium values 
are better than the charmonium ones. The trial function (80) seems somewhat better than (79), 
but the difference is not large, suggesting that the approximate variational solutions are reasonably 
accurate. Our results vary with our choice of parameter N, with N = 1.75 giving the best results 
(note that for our trial function (79), we must have > 1.5 to insure convergence of our integrals 
in (76); pushing down to this limit improves the results slightly. For (80), we need N > 0). 

C. NonAbelian terms present: (5 = 3,4, 5 cases 

For the (3 = 3, 4, and 5 cases, we follow the same procedure as in the previous section. Our 
equations become more complicated and one finds that it is easier to deal with the multidimensional 
calculations by resorting to Monte Carlo integration techniques [15]. Using the trial function (79), 
convergence is achieved much more quickly for the N = 3, 4 cases. Comparing these results with 
the (3 = 2 calculation for N = 3,4 is sufficient to give us an idea of the effect of increasing the 
number of gluons. From Table VII, for cc and N = 3, we can see that the (3 = 5 results are better 
than the (3 = 4 results, which are in turn an improvement over (3 = 3. However, none improve on 
the answers for (3 = 2. The same pattern holds for N = 4 and for bb. 
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VI. CONCLUSIONS 

Using the variational method and the Hamiltonian formalism of QCD, we have derived an infi- 
nite chain of equations for a quark-antiquark system interacting via an arbitrary number of gluons. 
These coupled equations are in principle exact, but not tractable. We attempt an approximation 
in which we decouple these equations using an ansatz borrowed from QED-type calculations. This 
leads to a sequence of (increasingly more complex) equations for the qq system with (3 = 0, 1,2, ... 
gluons present. NonAbelian effects appear only for (3 > 2. For (3 = 0, 1 we have only the Abelian 
Coulomb-type interaction present, and this leads to a Balmer-like mass-spectrum formula. We 
solve the (3 = 2 equation variationally, and work out predictions for the low-lying energy levels of 
charmonium and bottomonium (at least in the nonrelativistic limit). The results are encouraging 
since they show a substantial improvement over the Abelian approximation. We perform the same 
procedure for (3 = 3, 4, 5 but find that the results do not improve upon the ones for (3 = 2. It 
would be useful to see to what extent the results are improved if the relativistic versions of these 
calculations are performed. 
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APPENDIX A 
In Equation (31), we use the following definitions: 

Qi(fi,n, S2,s, A",k) = Qi(/3,n,s 2 ,s,0,k) + Qi(f3, n, s 2 , s, 3, k) 
Q 2 (/3,p,si,s,A",k) = Q 2 (/3,p,si,s,0,k) + Q 2 (/3,p,si,s,3,k) 

for A" = 0, and 

Qi(/?,n,s 2 ,s, A",k) = Qi(/3,n,s 2 ,s, A",k) 
Q 2 (/?,p,Si,s,A",k) = Q 2 (/3,p,si,s, A",k) 

for A" = 1,2. And, 

fe, j, X'", A', k') = X(0, k, j, 0, 0, k') + X((3, k, j, 3, 3, k' ) 

+ X(i3,k,j,0,3,k')+X(i3,k,j,3,0,k') for A',A / " = 
X((3,k,j,X"',X',k') = X(P,k,j,0,X',k') + X((3,k,j,3,X',k') for A' = 1,2; A' 
X(J3, k, j, X'", A', k') = X(J3, k, j, X 1 " , 0, k') + X(J3, k, j, X'" , 3, k') for A' = 0; A' 
X(P,k,j,X"',X',k')=X(P,k,j,X"',X',k') for A', A'" = 1,2 



where 



and 



QMn,s 2 ,s,X",k) = <^s 2 )^ (k) Hn-Ks) jf+1 



Q 2 (/3,p,si,s, A",k) 



\/w(n) oj(n — k) 

ii(p,fii) A"( k ) u(p-k, s) 
\/Mp) w(p - k) 



J, 



/3+1 



Q 3 (/3, fc,n,s 2 ,s) 



f(n,g 2 ) ^A fc (qfc) ^(n + qfc,g) 
yM n ) cj(n + q fc ) 



^3 



Q 4 (p, fc,p,si,s) = ====== J£ 

V w (p) ^(p + q/c) 



X((3,k,j,X'",X',k') = 2Kg(z k ,z j 



ev(k') • e Afc (q fc ) CAj(qj) ■ CA'"(qfc + q, - k 



CA fc (qjfc) • eA 3 (q,-) e A '( k ') 1 ey»(q k + - k') 



And, lastly, (summing over repeated indices), 
x(3\ {if 1 , tf 2 , •••if' 3 , } 

' <- tl *2 '2 »3 </3 ',9+1 ' 



Wl 1 



x(0 + l)l{fU Cis-C'J 
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J 2 +1 - {t Z i?i' 2 ■■■ tZ2 Q v } % i s h i s h+2 j j <% 



x((3 + iy.{t? lh t 



& p+ 

Zl 



p+1 ' 



J. 



• t Z2 if, 1 

l P 



t? , } 

'/3+1 l p+2 ' 



p+1 



0+1 



3 (**) = {*?V- 

x(/3-l)!{tf^ 2 ...^. 
x(/3-l)!{t- i2 .-.^. 

K$(z k ,Zj) = fazj d faz k e {t"P ^ • ■ ■ t* 2 t\} „ } ,v ; (5/^ 4 i( 5 I/3 + 1 



} tf v 5 hi 5 hj 5 Vii , Si 
*f / } 

t/3-1 1/5 J 
' } <*ili fysj fyj S l 



/3 + 1 

lp ' 



0+1 ' 



P "P+l ' 
Z P 



A closer look at our sums (A9) reveals the following relationships: 



/3+1 



J. 



/3+1 



^(^+1) = jI (^+0 

^ +1 = (/3 + l) jf +1 = (/3 + l) 2 Jf(^ +1 ) 
so that the expressions in (37) and (38) (i.e., 7(/3)) may be simplified to 



(A9) 



(AW) 



7(0) 



j/3+1 7/3 



<^4 ( z /3+l) "2 



J 



/3+1 7-/3 



J/3 7/3+1 



_ 9 ^f +1 ^4 ( Z /3+l) _ 9 



J/3 //3+1 
7/3+1 



( z /3+l) + - J 2 



Jo +1 J 4 (Zj3 + i) 



7/3 7/3+1 



7/3 7/3+1 



(/3 + 1) 3 ^' 



(All) 
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APPENDIX B 

Removing the permutation terms in (29) leads us to the set of equations: 



I P 4> f3 {p,n,s 1 ,s 2 ,q, A) 



uj p + uj n + m (mo — m) ( 1 ) + |Qfc| - E 

i / , f ,3, Qi(/3,n,s 2 ,s, A",k) „ +1 „ 
+ g m 2^2^^" / d fc /TO V (p,n-k,si,s,g,k,A,A 

A»=0 s J Vl k l 

2 ~ ^ 

i /- /" ^3? s i) s > A ,k) jij, „ 
# m /> /> Ca» / d A; = <p p+ (p-k,n,s,s 2 ,g,k,A,A 



A"=0 

/ QsiPi k, n, s 2 , s ) ,3-1/ , - r \ 

+ ? m L L 7r=^ v {P,n + q k ,s 1 ,s,q k ,X k ) 

s k=i vl^fel 

"9 m L L 7r=r ^ (p + q fc ,n,s,s 2 ,g fc ,A fe ) 

s fc=i Vl^fel 

-7 V ^((3^,3,1) <p p (p,n, Sl ,s 2 ,q k qj q h X k Xj Xi) — === r 

2 fe ^; Vlq*l I qui loil 

- -g- ^ fe, j, I, m) <f 4 (p, n, si, s 2 , q fc q., q ; q~ m , A,- A~ fe A; A m ) 

<5 3 (q fc +qj +q/ + q m ) 



I Qjl l<3Ul |qr. 



9 12 /• E 2 (P,k,j,l,\') (j) 13 2 (p,n,si,s 2 ,q fc + q ? - +q;,q fc q 3 q/,A',A~ fc A~- A/) 

■r- Ca' 2^ 



A'=0 fc/j^ 



V\^k\ |qul hi\ I qfc + qj + q«l 

d 3 ^ 



8 v,v*=o ^ Vlq fc Mq,Mk'l|qfc + q,-k'| 

x X(J3, k,j, X'", A', k') <^(p, n, ai, s 2 , q* + qu - k', k', q fc q~ , A'", A', A~ fc A 3 ) = (51) 



where the Q's and X are defined in APPENDIX A and 



for A", A' = 0, and 



f 2 (J3, k, j, X" ) = T 2 ((3,k, j, 0) + T 2 ((3, k, j, 3) 

F 2 (/?, fc, j, Z, A') = F 2 (0, k, j, 1, 0) + F 2 ((3, k, j, I, 3) (B2) 



f 2 (f3,k,j, X")=T 2 (f3,k,j,X") 

F 2 ((3,k,j,l,X') = F 2 ((3,k,j,l,X') (B3) 



for A", A' = 1,2. 
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Furthermore, 



Ti(P,k,j,l) = G% 3 (z k ,Zj,zi) e u Xk (q k ) e„ Ai (q ; ) Qfc • ^ (Oj) - |q*| £\A<\j) 



T 2 (P,k,j,X")=Gt 1 (z k ,z 



(q* + qj) A fc (qfc) [(qfc + qj) ■ e Xj (qj) - |q fc + qj | (qj) 
+ «A fc (Qfc) e ^ a, (qj) [q* • €a» (q* + Qj) Iqfcl e A" (q* + qj)] 
- «A fc (0* ) \" (q* + qj ) qfc • e A , (qj ) - |q* I e A , (Qj ) 



Fi(^,fc,j,Z,m) = {z^z^zuzm) e Afc „(q fe ) e^(qj) [e A ,(qj) e^ m (q m ) + e Ai (q/) -e Am (< 

F2(P,k,j,l,\') = K^~ 2 (z k ,Zj,zi) 

2 e Afc „(q*) e£ (q fe + q, + qj) [e^ (qj) (q,) + e A . (q,) • € A; (q,) 

+ 2 e A . „(qj) e^(q/) [e£,(q fe + q, + q/) e^ fc (q fc ) + £ A /(q fc + q, + q,) ■ e Afc (q fc )] 



Finally, (summing over repeated indices), 

x *fi * 2 *?» is " " " C Wi + Permutations [z u ■■■,z fS ] 
J i +1 = fZ i[i' 2 ■■■t* 2 *fjj *j i' *ii * ^3+2 j t' *i 



0+1 



# t Zl 



.2/3 



+ permutations[a, Z\,...,zp\ 



J 2 +1 ~ *ff • • • ^ *fjj J' *? i <*Ii t fys+a j J 



/3+1 



*?i «2 *?a i 8 ■ ■ ■ 2,3+2 + permutations [a, zi, . . . , zp\ 



.Zf3 



/3+1 



X 



h z k . . . + Z P 



t h-i h + P ermutations [ z i> ■ ■ ■ ,Zk, ■ ■ ■ ,zp] 



/3 + 1 



1/3-1 f/3 

x *fi i 2 ■ ■ ■ t Zk ■ ■ ■ 1? + permutations^, . . . ,z k , . . . ,zp] 
G%- 3 (z k , Zj , Zl )=i f Zk Zj Zl t\f v ^ • • • t Z2 tf/ 5 Vi j, 8 Vp+i i 8 h i S h _ 2 y 

x *fi ; 2 ' ' ' ^ 3 ' • • *^_ 3 «^_ 2 + permutationslzi , . . . ,z k ,Zj,z t , . 
G§-\z k , Zj )=i f Zk Zj c ■ ■ ■ i 22 ^ *ii i' i i fyu' 
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*fi h f h is'" tZk ^ " ' ^ + permutations[c, z u . . . ,z k ,Zj, . . . ,zp] 



-^2 { z k, Zj, Zl, Z m ) — fazi z k ft 



,, . , . t/ ■ ■ ■ t Z2 t z } <fy ,•/ ; Si' n _ i 8i x iSi 



tf 1 , ■■■ t Zk t Z i t Zl t Z 



t z /3 



+ permutations[zi, . . . ,z k ,Zj, zi,z m , ■ ■ ■ , zp\ 



K^~ 2 (z k , Z] ,zi) = f aZjZk faz ie tf iV2 ---t Z2 5 Vif S Vp+ii S hi S h _ ir 

tii h l hh ■■■ tZk tZi fZl ' ■■ t Z-2 h-i +permutations[e,z 1 ,. . . , z k , z"j, z u ■ ■ -,zp] 
K^(z k , Zj ) = f aZjd f aZke t^ Va ■■■t Z2 tf£ S l[f S^ +il 5 hi 5 h+lj/ 

*h * 2 *f a « 3 *?3 I* " " " *** ' " " Wi + Permutations [d, e, z x , . . . , z k z, , . . . , zp] 



(38) 
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Table I. Predicted energies for cc (nonAbelian terms suppressed) 



n 


(3 = 


(3=1 


(3 = 2 


Exp't(MeV) 


1 


a = 0.9542 


a = 2.181 


a = 2.672 


3096.88 


2 


m c = 1941 


m c = 1941 


m c = 1941 


3686 


3 


3795 


3795 


3795 


4040 


4 


3833 


3833 


3833 


4415 


5 


3851 


3851 


3851 




6 


3860 


3860 


3860 





Table II. Predicted energies for bb (nonAbelian terms suppressed) 



n 


(3 = 


(3 = 1 


(3 = 2 


Exp't(GeV) 


1 


a = 0.575135 


a = 1.31459 


a = 1.61038 


9.46037 


2 


m b = 5.10547 


m b = 5.10547 


m b = 5.10547 


10.02330 


3 


10.1275 


10.1275 


10.1275 


10.3553 


4 


10.1640 


10.1640 


10.1640 


10.5800 


5 


10.1809 


10.1809 


10.1809 


10.865 


6 


10.1901 


10.1901 


10.1901 


11.019 



Table III. Predicted energies for cc (nonAbelian terms present, (3 = 2) with trial function (79) 



n 


N = 1.75 


N = 2 


N = 3 


Exp't(MeV) 


1 


a = .7372 


a = 1.053 


a = 1.550 


3096.88 


2 


m c = 3370 


m c = 2717 


m c = 2267 


3686 


3 


3871 


3863 


3848 


4040 


4 


3961 


3948 


3923 


4415 


5 


4014 


3998 


3965 




6 


4050 


4031 


3993 




oo 


4223 


4190 


4121 





Table IV. Predicted energies for bb (nonAbelian terms present, (3 = 2) with trial function (79) 



n 


N = 1.75 


N = 2 


N = 3 


Exp't(GeV) 


1 


a = 0.5200 


a = .7016 


a = .9804 


9.46037 


2 


m b = 6.4710 


m b = 5.847 


m b = 5.4169 


10.02330 


3 


10.1999 


10.1928 


10.1781 


10.3553 


4 


10.2862 


10.2742 


10.2494 


10.5800 


5 


10.3372 


10.3219 


10.2902 


10.865 


6 


10.3710 


10.3533 


10.3166 


11.019 


oo 


10.5366 


10.5046 


10.4386 
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Table V. Predicted energies for cc (nonAbelian terms present, (3 = 2) with trial function (80) 



Tl 


N = .5 


N = 1 


N = 2 


Ern'fi MpV) 


1 
1 


m — 1 m A 


~, _ qi 07 


-V _ 1 KQO 

ui — i.jyz 


OUjU.OO 


2 


m c = 18119 


m c = 2934 


m c = 2244 


3686 


3 


3881 


3867 


3846 


4040 


4 


3979 


3954 


3920 


4415 


5 


4037 


4006 


3962 




6 


4076 


4040 


3989 




oo 


4270 


4205 


4114 





Table VI. Predicted energies for bb (nonAbelian terms present, (3 = 2) with trial function (80) 



n 


N = .5 


N = 1 


N = 2 


Exp't(GeV) 


1 


a = 0.0931 


a = .6258 


a = 1.0037 


9.46037 


2 


m b = 20.56 


m b = 6.0538 


m b = 5.3948 


10.02330 


3 


10.2098 


10.1960 


10.1767 


10.3553 


4 


10.3029 


10.2796 


10.2470 


10.5800 


5 


10.3586 


10.3288 


10.2871 


10.865 


6 


10.3958 


10.3612 


10.3130 


11.019 


oo 


10.5812 


10.5190 


10.4322 





Table VII. Predicted energies for cc (nonAbelian terms present, (3 = 2, 3,4, 5) with iV = 3 and 
trial function (79) 



n 


(3 = 5 


(3 = A 


(3 = 3 


(3 = 2 


Ex(MeV) 


1 


a = 1.886 - 


1.883 


a = 1.930 - 


1.928 


a = 1.799 - 


1.798 


a = 1.550 


3096.88 


2 


m c = 2179 - 


2180 


m c = 2155 - 


- 2156 


m c = 


2129 


m c = 2267 


3686 


3 


3841 - 


3842 




3839 




3836 


3848 


4040 


4 




3912 




3908 




3903 


3923 


4415 


5 




3951 




3946 




3940 


3965 




6 




3977 




3971 




3964 


3993 




oo 


4091 - 


4092 




4081 




4069 


4121 
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